* estimate quarter by quarter probits of being informed and plot fitted probs against actual status

quietly{
use WorkingData, clear

set more off

* check for non random assignment
xi: qui probit treatment  ladd i.dim i.area i.setton5 i.dp 
predict  PP, pr

* ************* GENERATE FIG A1 **************

collapse (p5) pp5=PP (p10) pp10=PP ( p25) pp25= PP  (p50) pp50= PP (p75) pp75= PP  (p90) pp90= PP (p95) pp95=PP , by(treat)
format pp* %3.2f
lab def lt 0 "Uninformed" 1 "Informed"
lab val treat lt
twoway rbar pp5 pp95 treat, barw(0.001) col(black)||rbar pp25 pp50 treat, barw(0.3) lcol(black) fcol(white) ||rbar pp50 pp75 treat, barw(0.3) lcol(black) fcol(white) xscale(range(-1 2)) xla(0 1, valuelab) legend(off) yla(, angle(0)) xti("Assignment status") yla(0.60(0.01)0.80, angle(0) labsi(*0.8)) graphregion(c(white))

graph export FigRndAssignment_A1.eps, replace

* ************* GENERATE FIG A1 : END **************

}

* ============================
* PRINT TO LOG DATA FOR FIG A1

list treat pp5 pp25 pp50 pp75 pp95
